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A system's invariance under Galilean transformation implies three locally conserved densities. 
Including them as variables, the thermodynamics is rendered explicitly frame independent, dissipa- 
tive mass currents are shown to vanish, and spontaneously broken Galilean symmetry becomes a 
sensible concept in condensed systems. 
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I. INTRODUCTION 

The invariance under Galilean boosts of a closed sys- 
tem implies the conservation of 

B = MRo, (1) 

where M is the system's mass, and Ro its center of mass 
coordinate at the time t = JlJ. Despite the appear- 
ance of an initial value that can always be set to zero by 
an appropriate choice of the coordinate, B is an additive 
and locally conserved quantity, in complete analogy to 
energy, mass, momentum and angular momentum. As 
these are all thermo- and hydrodynamic variables, it is 
a valid question why B is never included. This paper 
shows the usefulness, even necessity, of including B, sum- 
marized in the following list. 

i) The inclusion of B as a variable is a necessary condi- 
tion for the formulation of a frame independent thermo- 
dynamic theory: Starting from the rest frame expression, 
say dE — TdS, it may seem all right to add V- dG for a 
boosted system with the momentum G, and add Q ■ e?L 
if it also rotates with the angular momentum L, but this 
is not enough. As shown in section this is only correct 
for frames in which f2 1| G, while generally the energy E 
also depends on B. 

Before going on with the list, it is convenient to intro- 
duce a name for 

B = M(R-Rt) = Jd 3 x(gx-gt), (2) 

where R = Ro + Rt is the time dependent center of mass 
coordinate, with a constant time derivative R = G/M, 
while g and g denote the density of mass and momentum, 
respectively. Following Schwinger Q , we shall refer to B 
as the "booster" , and to b = gx—g t as the "booster den- 
sity" . Including the energy's dependence on the booster, 
dE = . . . + A- dH, the conjugate variable A will turn out 
to be ft x R. 

ii) One of the more direct results of the thermody- 
namic theory are the equilibrium conditions. Again, the 
correct derivation requires the inclusion of the booster, 
although this is never done. If one derives them follow- 
ing Landau and Lifshitz || by maximizing the entropy 
while holding constant: the energy, mass, momentum, 



and angular momentum, we find as equilibrium condi- 
tions: constant temperature and chemical potential, with 
the velocity given asv = \ + flxx, where V and fl 
are again constant. Due to the lack of a time dependent 
term, however, this velocity is less general than the usual 
expression for a rigid body, v — R + i~i x (x — R), or 
equivalently 

v = (R - ft x R ) + n x x - (SI x R) t . (3) 

(The difference is in the last term.) This is worrisome as 
there is no reason whatever why a system executing the 
general motion should not be in equilibrium. If, however, 
the booster is also held constant in the above calculation, 
one does arrive at the general expression, Eq(||). At the 
same time, one finds that the chemical potential, instead 
of being constant, now satisfies dtv + V/i = 0. Reas- 
suringly, this is exactly the expression the Navier-Stokes 
equation reduces to for vanishing entropy production, in 
equilibrium. 

iii) Local conservation of the mass density holds if the 
continuity equation, dtQ + V- j = 0, is satisfied — irre- 
spective what form the mass current j actually assumes. 
Usually, this form is taken as j — gv, although it has 
never been properly deduced; rather, it is accepted as 
a statement of microscopic plausibility, or the summary 
of countless experiments. So no objection was, or could 
have been, raised, when Dzyaloshinskii and Volovik, in 
their classic paper , propose to include dissipative mass 
currents such as 

j-gv^V/j,. (4) 

(These dissipative terms results from a tempting, even 
natural, step to take when setting up the hydrodynamic 
equations.) On the other hand, there is a well-hidden 
footnote by Landau and Lifshitz J5| that purports to rule 
out this type of terms — but actually falls short of being 
ironclad if scrutinized: It invokes the center of mass mo- 
tion to show that Jd 3 xj = Jd 3 xg must prevail, where 
the integration is to be taken over the volume of the 
system. For reasons of Galilean invariance, the momen- 
tum density is given as g = gv. So a clear-cut proof 
of j = g would indeed serve as a sound argument for 
ruling out any dissipative mass currents. Unfortunately, 
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Jd 3 xj = Jd 3 xg is less confining, and the reader is left 
wondering about terms that vanish only if integrated. 

As will be discussed in section [II, taking the booster- 
density q x — g t as a locally conserved quantity that sat- 
isfies a continuity equation quickly leads to the result 
j = g. The proof takes place in very much the same 
way as that deducing the symmetry of the stress tensor 
from the local conservation of angular momentum. (This 
result also has relativistic ramifications, as it shows that 
out of the two versions of relativistic hydrodynamics, by 
Eckart and by Landau-Lifshitz, only the latter is a proper 
generalization || . Any linear combination of both is also 
ruled out.) 

iv) The angular momentum frequently contains an ex- 
tensive part, s, that is usually referred to as its intrinsic 
contribution, 



L = Jd 3 x(xxg + s). 



(5) 



One example of s is the spin density. Due to the close 
relationship between L and B — the latter being, rela- 
tivistically, part of the 4-tensor of the angular momen- 
tum density — it is certainly not far-fetched to question 
whether a corresponding extensive contribution to the 
booster may exist, 



B = jd 3 x(gx - gt + k). 



(6) 



(Relativistically, the components of the 4-angular- 
momentum density mix under a Lorentz boost. And s, 
the extensive contribution of the angular momentum in 
one frame, will lead to k in others. The more precise 
question is then: Whether k exists in the rest frame, 
same as s.) Microscopically, k may be conceived as aris- 
ing from individual mass dipole moments of the molecules 
or atoms that are not compensated by shifting the coor- 
dinates of the particles, say because they are embedded 
in a lattice, pivoted off center. 

Although the intrinsic angular momentum s has fre- 
quently been included in hydrodynamic considerations, 
other authors see this as ad hoc and inconsistent, be- 
cause s is a relaxing variable, and its inclusion only affects 
high frequency, non-hydrodynamic phenomena. Redefin- 
ing the momentum density, they then argue that s can 
always be chosen to vanish, and therefore never needs 
to be considered j7| . While the first part of this stance 
is fairly convincing, the final conclusion excludes a whole 
category of hydrodynamic phenomena. These include es- 
pecially the Einstein-de Haas effect, in which a station- 
ary magnetizable body starts to rotate when the external 
magnetic field is turned off [|| : The vanishing velocity v 
of a stationary body also compels the momentum density 
g = gv and the orbital angular momentum density x x g 
to be zero. But the total angular momentum L may still 
be finite in the presence of a magnetic field H if there is a 
finite intrinsic angular momentum s ~ 7i. If Ti. is turned 
off, s vanishes, yet L must remain constant. Hence the 
system starts rotating, as observed, to compensate. 



The analogous effect for the booster would be given 
by aligning the microscopic mass dipoles with an exter- 
nal electric field, ie k ~ S. Turning the field off gets 
rid of the alignment and kills k. The orbital part MRo 
compensates by displacing the crystal. 

The proper way to account for this type of effects lies 
in the middle ground between the above two extreme 
points of view — taking s as either completely indepen- 
dent or utterly negligible. As discussed in section |y|, we 
should in fact accept s as an thermodynamic variable, but 
not a hydrodynamic one (the precise meaning of which 
will become clear there). The close analogy between s 
and k makes it easy to treat the latter along the same 
line. One of the results, based on fairly general thermo- 
dynamic considerations, is that the booster's analogue of 
the Einstein-de Haas effect does not exist. 

v) Finally, a truly serendipitous result: the collective, 
hydrodynamic behavior of systems that spontaneously 
break Galilean symmetries. This subject has until now 
eluded clarification because the relevant conserved quan- 
tity, the booster, has been neglected. Again, the analogy 
between the angular momentum and the booster comes 
in handy here, as broken rotational symmetry is a well 
understood concept |t|J|. The order parameter u of the 
broken Galilean symmetry obeys the equation of motion, 
dtVL + V/i = 0. Despite the obvious similarity to the 
Josephson equation [jl0| , the Goldstone modes of broken 
Galilean symmetries are not propagating second sound 
modes, but resemble the orbital waves of nematic liquid 
crystals. 

Unfortunately, we have nothing to say concerning 
possible microscopic mechanisms that realize broken 
Galilean symmetry. The system that spontaneously 
breaks this symmetry needs to display indifference with 
respect to the inertial system it chooses while being able 
to sustain a velocity gradient in equilibrium. And the 
relevant order parameter u is a velocity difference. Such 
as system is admittedly hard to conceive. 



II. FRAME-INDEPENDENT 
THERMODYNAMICS 

Let us begin by deducing the thermodynamic equi- 
librium conditions. Maximizing the entropy, or equiv- 
alently, minimizing the energy, while holding constant 
the entropy and the conserved quantities, including es- 
pecially the booster B, we have 

S[jd 3 xe~f Jd 3 xs~fi fd 3 xg-V-fd 3 xg 

-Cl-Jd 3 x{xxg) - A-fd 3 x{gx -gt)] = 0, (7) 

with T, (m, V, f2, and A being constant Lagrange para- 
meters. Employing the local Gibbs relation, 



de = T ds + /i dg + v ■ dg . 
this can be written as 



(8) 
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/ d 3 x [ (T - f) Ss + (ji - fl - A-x) Sg 

+(v- V -£lxx + At)-8g] =0. (9) 

From this follow the Euler-Lagrange equations 

T = f, a = fi + A-x, v = V + (lxx-At. (10) 

Comparing the last of Eqs(10) to the velocity field of 
Eq(H), we find complete agreement with 



A = f2xR, n = £l, V = R fixR . 



(11) 



Excluding the conservation of the booster is equivalent 
to setting A = £1 x R = 0, cf Eq(@). This leads to the 
results of Landau and Lifshitz, mentioned in the intro- 
duction, correct for frames in which Q || R. 

In a more deliberate approach, the input of Eq(|3|) is 
not necessary: Varying also the time and space coordi- 
nates while minimizing the energy yields two additional 
terms, 



fd 3 x 



(nxg-Ag)-6x + (A-g)5t\=0. (12) 



Taking uniform transformations in time and space, 
St, 5x = const, as independent variations, we conclude 
A = fl x R and (redundantly) A R = 0. 
Locally expressed, Eqs(|l(]) become 



VT = 0, 



d tVj + djVi = 0, d t v + Vu = 0, (13) 



which are useful for the off-equilibrium considerations be- 
low. There are three corollary conclusions to draw: 

(i) Integrating Eq(|8j) over the system's constant volume 
while heeding Eqs(M), we obtain the Galilean covariant 
version of the extensive, basic formula of the thermody- 



dE = fdS + ftdM + V- dG + Ct ■ dL + A- dB. (14) 

If the volume is allowed to vary, there is an additional 
term, — §d 2 x P du, with P = ag + Ts + v ■ g — e and du 
denoting the displacement of the surface, along the sur- 
face normal and at the area element d 2 x. The unusual 
form is related to the fact that P is not a constant, and 
the energy change depends on where the volume change 
takes place. For P — const, this term reduces to the 
usual form —P§d 2 xdu = —PdV. The last three terms 



of Eq(14), ie the kinetic part of the energy, can be inte- 
grated (for a sphere) to become 



G 2 /2M + (L - R x G) 2 /26, 



(15) 



as it should. O is the moment of inertia in the center of 
mass frame. 

(ii) Despite the manifest frame dependence of \x and 
v in Eq([l0|), most quantities are properly frame inde- 
pendent. Take eg the density distribution under rota- 
tion: The chemical potential of the local rest frame, 
fj,o = [i + \v 2 , is a function of T and g, hence V/io = 



(du Q /dg)'Vg for constant temperatures. On the other 
hand, we deduce from EqsdlG, Hi]) that 



V/iQ = V(/i 



\v 2 ) 



±V[f2 x (z-R)] 2 



(16) 



depends only on the velocity in the center of mass frame, 
and not on R, the center of mass velocity of a given 
frame. So V/^o, and therefore Vg, remain unchanged 
under a Galilean boost. 

(iii) The Navier-Stokes equation is dv/dt + VP/ g = 
for diVj + djVi = 0, and dv/dt + V/xo = if in addition 
VT = 0. Rewriting dv/dt = d t v + (v-V)v = d t v-\Vv 2 , 
it finally reduces to the third of Eqs(|13|), d t v + V fi = 0. 



III. THE DISSIPATIVE MASS CURRENT 

The inclusion of the booster is clearly important for 
static, thermodynamic considerations, but it is equally 
relevant for off-equilibrium, dynamic situations: First, 
the proof for j = g, or that no dissipative mass current 
is allowed. As this proof follows closely the one that de- 
duces the symmetry of the stress tensor from local con- 
servation of angular momentum, it aids comprehension 
to present both simultaneously. Rewriting the continuity 
equations for mass and momentum, 



as 



d t g + diji = 0, d t gi + djllij = 0, 



d t bi + dj(xijj - Uijt) = ji - gt 



(17) 



(18) 
(19) 



where £ = xxg and b = gx — gt, we immediately see 
that eijfellfcj and j — g have to vanish for the angular 
momentum and the booster to be locally conserved. 

Now, one may argue that it would be quite enough if 
these two expressions can be written as divergences of 
some currents, e ijk Tlkj = d k a ik and j 4 - g t = d k J ik , so 
they need not vanish. But this neglects the following two 
points: First, qualitatively, I and b contain the reference 
to the origin of the coordinate, g and g do not. Con- 
sider a small volume element far away from the origin, 
at distance R and t = 0, then £ and b = gx will scale 
with R. Now, as we can change R by simply relocating 
the origin, without altering any of the physics, d t £, d t h 
and their respective fluxes must also scale with R. On 
the other hand, because Hij and j (as fluxes of g and g) 
do not, neither do e^kHkj = d k a ik and ji - cji = d k J ik , 
which therefore must vanish. 

Second, the quantitative argument. Being currents of 
hydrodynamic variables, and ji are themselves func- 
tions of these variables, and of their spatial derivatives. 
To lowest order, with no spatial derivatives, Iljj and ji 
cannot possibly be written as divergences of some cur- 
rents. In the next order, the terms are mostly dissipa- 
tive. Taking eij k H k j, j — g as finite, we find the entropy 
production to be given as @Jll]] 
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K 



e ijk IL jk Sli-(j-g)-VtJ,, ft = ±Vxu. (20) 



According to the rules of irreversible thermodynamics, 
if the entropy production is given as a sum of products, 
1Z = TiXiYi, we may take the fluxes Xi as proportional 
to the thermodynamic forces Yi |yj, hence 



e-yfelfj-fe — C(i) Qi, 3 - 



"C(2) V/z, 



(21) 



where and £ (2) are transport coefficients, similar to 
the viscosity. (More accurately, the Xi are a linear com- 
bination of the Yi, so off-diagonal, cross terms are also 
possible. However, these terms always vanish if the di- 
agonal ones do, so they need no extra consideration.) 
Eqs(pT|) lead to the entropy production 



K 



(22) 



which, however, contradicts Eqs(|l3|): ft and V/i need 
not vanish even if all the three fields there do, in which 
case equilibrium reigns with 1Z = 0. Therefore £ (2) 
are zero, and e^Hy, J — g vanish. 



IV. INTRINSIC BOOSTER-DENSITY 

Allowing now intrinsic contributions, as in Eqs(|, 
the energy density depends on two additional variab 
s and k, or 



es, 



de = Tds + n dg + v ■ dg + u> ■ ds + a ■ dk . 



(23) 



Barring an instability, the energy e is minimal for s, k = 
0. Therefore, an expansion around the minimal value, 
Smini would usually yield a quadratic dependence, e = 
£min + ^(l(i) s2 + 7<2)k 2 ). (Isotropy is assumed, otherwise 
7d)i 7(2) are tensors rather than scalars.) Being deriva- 
tives, oj and a assume the form 



w = 7(1,8, a = 7 (2) k. 



(24) 



If the magnetic field is included, de = . . . + Ti ■ dB, the 
energy is minimal for s, k, B — 0, and the same expansion 
(again for an isotropic medium) yields 



w = 7 ( i,s + 7 C 3)B. 



(25) 



It is this cross dependency that leads to the Einstein-de 
Haas effect discussed in the introduction ||. Similarly, 
an electric field leads to 



a = 7 (2) k + 7 (4) X>. 



(26) 



Because t + s and b + k are locally conserved quantities 
and satisfy continuity equations, the new variables obey 
the equations of motion, 

d t Si + djXij = ei jk U jk , dth + dj^ij = g l - j t , (27) 

where and 7^ are the respective currents, functions 
of the local hydrodynamic variables, the explicit form of 



which we shall not derive here. The source terms on 
the right hand side cancel those of Eqs([l9|), and ensure 
that £ + s and b + k satisfy continuity equations. (The 
vanishing of ejjfcHjfc and g — j was deduced in the last 
section as a direct consequence of t and b being locally 
conserved quantities. As this is no longer true, we have 
no prior knowledge of whether tijk^ijk and g — j are fi- 
nite.) Going through the same hydrodynamic procedure 
as that leading to Eq(p0|), we find the modified entropy 
production, 



n 



+ ey fc II 3 -*(fti - uJi) -{3- ff)-(V/i - a), (28) 



and conclude as before, 

eijkRjk = C<i)(^i - Wj), 3 - g = -C ( 2)(V/i - a), (29) 

If the two coefficients were exceedingly large, , C<2> - * 
00, we may divide each of Eqs(^) by the respective co- 
efficient and find the equations substituted by the two 
constrains, 



uj = £1 , a = V/i . 



(30) 



The first implies instantaneous establishment of partial 
equilibrium by local exchanges of angular momentum, 
between the reservoirs £ and s; the second implies par- 
tial equilibrium by local exchanges of booster, between 
b and k. So instead of an independent dynamics, s and 
k assume the appropriate values satisfying Eqs(|30|), in- 
stantaneously during hydrodynamic processes. 

For generic circumstances, this is indeed the case. 
There are two terms in Eqs(p9|), 

C ( i)W = C(i)7(i)S = s/r (1) , C ( 2>a = C( 2 )7( 2) k = k/r (2) , 

which show that both Eqs(^) are relaxation equations, 
and that Eqs(|30|) hold after the time t {1) and r (2) , respec- 
tively. Generally speaking, these characteristic times, af- 
ter which local equilibrium is established, are microscopic 
in nature and much faster than typical hydrodynamic 
time scales. And taking these times as negligibly small 
is equivalent to the limit £ (2) — > 00, and equivalent 
to substituting Eqs(|30|) for the equations of motion (p7|). 
Note, however, that this dynamic dependency does not 
at all mean s and k can not remain full-fledged ther- 
modynamic variables, as Eqs( p3| , p4 25 2^) clearly remain 
valid. 

From a more elevated point of view, we understand 
that the universality and simplicity of hydrodynamic the- 
ories are achieved by eliminating all relaxing variables, 
and expunging all relaxation equations. As a result, the 
hydrodynamic theory is confined to frequencies much less 
than all relaxation rates 1/T (i) . 

(A possible source of confusion is the seemingly odd 
fact that the transport coefficients Q 1} and C<2> were de- 
duced to be zero in the last section, yet argued to be 
diverging here. The explanation lies in the difference 
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between the calibers of the arguments: In the last sec- 
tion, the coefficients were not negligibly small. Rather, 
they had to vanish identically to not contradict thermo- 
dynamics. Here, the coefficients are finite, but well ap- 
proximated by oo when compared to hydrodynamic time 
scales. In other words, if they are finite, they are very 
large - and there is no contradiction.) 

One problem remains: The two terms eijk^-jk and j — g 
are now indeterminate, as they are both products of an 
infinite and a vanishing factor, cf Eqs(|29|). Nevertheless, 
they need to be known before we may make use of the 
continuity equations for mass and momentum, Eqs([l7|). 
A clever way to accomplish this is to go over to the new 
densities 0, 



<7=9 + 5( Vxs )' £> = £-V-k. 



(31) 



They are obtained by starting from 
incorporated, 



E3), with Eqs(P 



de = . . . v ■ dg + Q • ds + /i dg + V/i • dk. (32) 

This is partially integrated to yield 

de = . . . + v ■ dg + fidg, (33) 

which shows first of all that the energy density is well 
accounted for by the new densities. Putting the surface 
slightly beyond the system's volume, we also have 



Jd 3 x g = J d 3 x g, Jd 3 x g = J d 3 x g 



(34) 



as the respective second term of Eq(|31|) are surface terms. 
So we conclude that g and g are also valid densities for 
momentum and mass. Most crucially, because of 



and 



/ d 3 x (x X g) = Jd 3 x (x x g + s), (35) 



fd 3 x(gx -gt) = fd 3 x{gx -gt + k), (36) 



these new densities clearly absorb the intrinsic contribu- 
tions, s, k, and therefore have fluxes for which 



eijkKjk =0, J 



0, 



(37) 



hold. So it is indeed a good idea to take g and g as the 
hydrodynamic variables. (But it is now also possible to 
return to the original variables: Given Ily and J, it is 
easy to obtain Ily = fLy + \e i3 kd t s k , 3 = 3 <9*k.) 

All this does not mean that no experimental conse- 
quences are related to s and k, as claimed in j?]]. For 
instance, the constitutive relations, g = gv and Eq(25), 



g = gv + ^Vx s(n,B), 



(38) 



can not be without consequences in the Navier-Stokes 
equation, dtgi = TLij: A temporally varying but spa- 
tially homogeneous B-field will lead to an acceleration, 



dt{&v) = -9 t [Vx s(B)/2], at the system's outer rim, 
the only place where V x s is nonzero. The resultant 
rotary motion there then propagates into the bulk, elas- 
tically in a solid and viscously in a liquid. Clearly, this is 
a temporally resolved, hydrodynamic description of the 
Einstein-de Haas effect. 

Other thermodynamic cross dependencies of s will lead 
to similar behavior. The diagonal term all by itself, how- 
ever, while necessary for maintaining the thermodynamic 
stability, is probably always negligibly small, as 

Vx s/2 = Vx Sl/2j m 

is of the same order as the viscous terms. 
In spite of the one-to-one analogy, 

d t Q = d t [g- V-k] 

instead of Eq(|38|), there is one crucial difference between 
s and k: The quantity 

de/dk = V/u = O x R 

depends on the inertial frame. As a result, thermody- 
namic cross dependencies are ruled out for k. To un- 
derstand why, consider k(£), the electric analog of s(B). 
Due to the Maxwell relation. 



dkj/dSi = dVi/d{djn), 



(39) 



that follows from the energy de — . . .SidT>i + {djji)dkj, 
k's dependence on £ implies X>'s on V/i, a quantity that 
can be made to vanish by choosing an appropriate frame. 
In the absence of a magnetic field, T> is essentially frame 
independent, X>(V/x) = X>(0), or dVjd^d^i) = 0. Hence 
dkj/BSt = 0. 

The diagonal term 7 (2 >k = V/i = f2 x R remains, 
but the associated hydrodynamic effects, <9tV • k ~ 
(dn/dg)dtV 2 g are again an order smaller than the dis- 
sipative terms. They may probably always be ignored. 
Nevertheless, understanding the treatment of a finite k 
is useful, both for its own sake and for the considera- 
tion below of broken Galilean symmetry, (although these 
system do not necessarily have a finite k.) 



V. SPONTANEOUSLY BROKEN GALILEAN 
SYMMETRY 

Thermo- and hydrodynamic theories account for the 
generic behavior of macroscopic systems that results from 
conservation laws and broken symmetries. Although the 
specificity of its predictive power is, compared to the gen- 
eral input, frequently amazing, some subjects are simply 
beyond its reach. One example is the question concerning 
possible microscopic mechanisms which lead to broken 
Galilean symmetries, or in fact whether it at all exists. 
Fortunately, we do not need this knowledge to give a full 
account of the collective behavior of a system that breaks 
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this symmetry, and hereby filling the macroscopic half of 
a large gap in one of our more basic concepts. 

Starting with a system that breaks rotational invari- 
ance in all three directions, (say biaxial nematics or the 
B-phase of superfluid 3 He uM,) we take the infinitesi- 
mal rotation angle d9i as the order parameter, with the 
energy depending on its gradient, 



de = 



Qidsi + 4>ijd(dj6i). 



(40) 



An expansion of the energy e leads to <pki ~ djOi. The 
equations of motion in the rest frame possess the usual 
Hamiltonian form, 



d t Ui + (vkdk)ui + difi - CdjJij = 0, 
d t g + digi = 0, g. t = gvi- djjy. 



(44) 
(45) 



In conjunction with the rest of the hydrodynamic equa- 
tions that essentially retain their form from the isotropic 
liquid, we find the sound mode to be unchanged, while 
three additional modes, of coupled temperature and u- 
motion, are (as in the orbital case) of the form lu ~ q 2 , 
with the coefficients being complex. 



d t 6i = Se/Ssi = Hi, d t Si = -Se/S8i = dj 



(41) 



Eliminating dtSi as outlined in the last section, the con- 
tinuity equation for the momentum gets modified: 



dtfji = dt9i 



1 



djdfSk — -£ijkdjd m cj)km 



which leads to three pairs of orbital waves, of the form 
u> ~ q 2 , with complex coefficients |0|. (A slight rear- 
rangement is needed to render the seemingly antisym- 
metric stress tensor symmetric.) 

The order parameter of broken Galilean symmetry is 
a velocity field, Uj, with the energy depending on its gra- 
dient, 



de = 



(di(i)dfa + Jijd^djUi). 



(42) 



An expansion leads to Jki ~ djUi. The equations of mo- 
tion in the rest frame again possess the usual Hamiltonian 
form, 

d t Ui = -8e/5ki = — dth = Se/Sui = -d 3 ;J i5 ■. (43) 

Eliminating dtki as outlined in the last section, the con- 
tinuity equation for the mass density gets modified: 

d t Q = d t Q - didth = —di(gvi - djJij), 

where the total mass current equals the momentum den- 
sity. Including nonlinear terms and the diagonal dissipa- 
tive one, we have (dropping "), 
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